Nonadibaticity and single electron transport driven by surface acoustic waves 
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Single-electron transport driven by surface acoustic waves (SAW) through a narrow constriction, 
formed in two-dimensional electron gas, is studied theoretically. Due to long-range Coulomb interac- 
tion, the tunneling coupling between the electron gas and the moving minimum of the SAW-induced 
potential rapidly decays with time. As a result, nonadiabaticiy sets a limit for the accuracy of the 
quantization of acoustoelectric current. 
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Recently, a new type of single electron devices was in- 
troduced. In the experimentstJ, surface acoustic waves 
(SAW) induce, via piezo-electric coupling, charge trans- 
port through a point contact in GaAs heterostructure. 
When the point contact is biased beyond the pinch-ofF, 
the acoustoelectric current develops plateaus, where 



/ = N^ef. 



(1) 



Here / is SAW frequency, and iVg is an integer. The 
plateaus were shown to be stable over a range of tem- 
perature, gate voltages, SAW power, and source-drain 
voltages. Remarkably high accuracy of the quantization 
( l|) , and high frequency of operation {f 3 GHz) imme- 
diately suggest a possibility of metrological applications 
of the effectfl. 

However, deep understanding of these results is still 
lacking. Qualitatively, the effect is explained by a sim- 
ple picture of moving quantum dotscl. Electrons, cap- 
tured in the local potential minima ('dots'), created by 
SAW, are dragged through the potential barrier. The 
strong Coulomb repulsion prevents excess occupation of 
the dot. Increase of the SAW power deepens the dots, 
more states become available for the electrons to occupy, 
and new plateaus appear. By changing the gate voltage, 
the slope of the potential barrier can be lowered, which 
has a similar effect. 

Interestingly, the quantization was not observed in the 
open channel regimeB, although it shoiild be expected 
on the quite general theoretical ground[_ For the mech- 
anism of the quantization, discussed inci, it is essential 
that the DC conductance for each instantaneous configu- 
ration of the SAW-induced potential is zero. In the open 
channel regime it would require the channel length to be 
much longer than SAW wavelength A, which is difhcult 
to realize. In the experimentstJ this problem is avoided, 
since in the pinch-off regime the DC conductance is zero. 
However, as explained below, the rapid change of SAW 
potential near the entrance to the channel creates a new 
trouble, leading to nonadiabatic corrections to d|). 
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FIG. 1. Geometry of the system: (a) split gates; (b) re- 
gions, occupied by two-dimensional electron gas; (c) depleted 
region. The arrow indicates direction of SAW propagation. 



Long-range Coulomb interaction plays crucial role in 
this phenomenon. The two-dimensional electron gas 
(2DEG) is depleted in the vicinity of the gates (see 
Fig.l). On the other hand, in the depleted region screen- 
ing is lacking. The important parameters of the prob- 
lem can be understood from the solution of the electro- 
static model. In this approach, one assumes that, since 
the Fermi velocity uj? is large compared to the sound 
velocity Vg, the 2DEG is able to follow the changing 
in time SAW-induced potential. Therefore, it is sufh- 
cient to consider an instantaneous electrostatic problem, 
treating time as a parameter. Since the screening length 
in 2DEG (~ 10 nm) is much smaller, than A ~ lO'^ nm, 
one can assume that the SAW-induced potential is com- 
pletely screened in the 2DEG-occupied region. There- 
fore, one has to solve the Poisson equation subjected to 
complicated boundary conditions. In particular, the po- 
tential at the gates {(p = Vg), as well as the potential 
of 2DEG region {ip — 0), and density p = pg of the 
positive background charge in the depleted region are 
fixed. Furthermore, for simplicity, one can take an effect 
of SAW into account through a weak periodic modula- 
tion of p: p — > po + {x, t). The self-consistent solution 
would yield the location of the edge of 2DEG, potential 
in the depleted region ip (r) , and number density n (r) 
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of 2DEG. Then, using the Thomas- Fermi- type relation 
U (r) -I- (ttTt?' / m*^ n (r) = ep, one would be able to find 
an effective confining potential U (r). However, the de- 
tails of the full solution are not required, since the most 
important properties £an be understood from the follow- 
ing simple arguments^. 

Firstly, it is clear that the potential ip (a;, y) has a min- 
imum in the gap between the gates (see Fig.l). Secondly, 
the very presence of plateaus shows that the charge states 
of the dot are separated by the finite energy gaps. Since 
both A^o = odd and A'o = even plateaus are observed the 
energy gap is associated with Coulomb repulsion rather 
than single-particle level spacing. This means that for 
any given x, ip{x,y) has a sharp minimum near y = 0. 
The shape of this minimum is slowly changing with x, 
and this change is controlled by the geometry of the de- 
vice. Note that the smoothness of the change of the con- 
fining potential in y-direction is supported by the fact 
that the same systems exhibit very nice pattern o£ con- 
ductance quantization in the open channel regimecl. Fi- 
nally, the weak perturbation of the background charge 
density \Sp\ <^ po should not affect significantly the po- 
sition Xo of the 2DEG edge. 
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FIG. 2. The potential along the line y = 0. The dashed 
line shows the effective potential, that confines 2DEG to the 
region x < xo- The SAW-induced potential minimum (a;i) 
moves to the right with the sound velocity Vs- 



As function of time, Xq (t) oscillates with SAW fre- 
quency /. However, the amplitude of the oscillations is 
small, compared to A, therefore the velocity of the edge 
is negligible compared to Vs- On the other hand, the 
SAW-induced potential minimum (the dot) moves away 
from the edge with precisely the sound velocity (see 
Fig. 2). Therefore, the width of the potential barrier, sep- 
arating 2DEG and the dot, increases linearly with time. 
This in turn means that the tunneling coupling between 
2DEG and the level, localized in the dot, rapidly de- 
creases, approximately exponentially. The characteristic 
time can be estimated as t ~ Iq/vs, where is the dis- 
tance over which the localized wave extends under the 
barrier. Since, evidently, /q ^ A, the relation 



fT « 1 (2) 

holds. Other parameters, such as the height of the bar- 
rier W, the energy of the locaHzed level eo, etc., change 
during the time, which is of the order of SAW period 
1//. Thus, the time dependence of all these parameters 
during the time r can be neglected. 

Due to the rapid decrease of the tunneling coupling, 
the thermal equilibrium in the system can not be main- 
tained, causing fluctuations of the occupation number of 
the dot. This results in nonadiabatic corrections to the 
quantized values of the acoustoelectric current. Whether 
these corrections have a significant impact on the accu- 
racy of the quantization, depends on the value of the 
characteristic energy scale Ti/t^ as compared to other en- 
ergy scales in the problem. To obtain an order of mag- 
nitude estimate, we expand the potential near the mini- 
mum xi (see Fig.2), V {x) ~ Aq^ {x - xi)^ /2, q = 27r/A. 
The amplitude A is related to the single particle level 
spacing in the dot A via Aq^ = m* (A/Ti)'^ (m* is effec- 
tive electron mass), to the 'size' r of the localized wave 
function via Aq^r^ ^ A, and to the charging energy Ec 
via Ec ~ e^/er. Iq is estimated from WKB relation 
[Tl/Iq) ~ 2m*W. Assuming that W ^ A, it gives us 
four equations for five unknown quantities. Additional 
relation follows from the experimental results. It was 
demonstratedu, that the quantization disappears above 
the activation temperature T* ^ 10 K, which we identify 
with the charging energy. Using typical parameters for 
the experimentsEI, we find t ^ lOpaj. All the param- 
eters manage to pass minimal consistency requirements 
r ^ A, A <^ A, Ec, fr <^1. Since the corresponding en- 
ergy scale Ti/t ^ 0.1 meV, the nonadiabatic effects may 
have significant influence on (^ at low temperature. 

This can be understood from the following model 
Hamiltonian: 

H — Heg + Hdot + Ht- (3) 

Here 

^<=9 = H ^fe^aCfe^ (4) 
ka 

describes electron gas in the lead, 

Hdot = Endi^dna + Ec (N - Ngf (5) 

is the Hamiltonian of the dot [N — X^no- '^no-'^no- is the 
total number of electrons in the dot), and 

HT = V{t)Y,Aadna+'R-C., (6) 

fcncr 

describes the tunneling coupling with time-dependent 
tunneling amplitude. We have included only one lead in 
the model, since tunneling coupling to the second lead is 
negligible. The electron gas in the lead is assumed to be 
in thermodynamic equilibrium at all times, by virtue of 
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the inequahty vp ^ Vs- As discussed above, we neglected 
time dependence of various parameters in (|5|), due to the 
separation of the time scales /r <C 1. The last term in 
(g) describes the intra-dot Coulomb interaction. The pa- 
rameter Mg describes the effect of the gate voltage. Since 
the width of the plateaus is approximately independent 
on the plateau's number, it is a good approximation to 
assume that A/" is a linear function of Vg . The most im- 
portant ingredient of (|^) is the time-dependent tunneling 
amplitude. We take it in the form 



(7) 



(other possible choices will be discussed below). The 
time-independent version of (|^^ is commonly used in 
the theory of the Coulomb blockade. Similar models have 
been also employed to study transfer of charge during 
the atom-surface. scatteringQ, and nonadiabatic effects in 
charge pumpingJ. 

Given that the system, described by (|[^), is in thermo- 
dynamic equilibrium at t — — oo, our task is to calculate 
the occupation of the dot at t — > oo, Nq = {N)^^ The 



acoustoelectric current is related to Nq through (y 

Away from the Coulomb blockade degeneracy points 
(half-integer Ag), when the inequality 



2Ec lA'g - no - 1/2| > max{r, 1/t} 



(8) 



is satisfied, the time-dependence is too slow to cause the 
transitions between different charge states of the dot. In 
(^, no is integer part of Afg] units where h = ks = I are 
used throughout the rest of the paper. In this respect 
the evolution of the system is almost adiabatic, and the 
occupation of the dot A^o is expected to coincide with 
equilibrium occupation, corresponding to the Hamilto- 
nian (^), with T replaced by the effective temperature 
Teff ~ max{T, l/r}. However, in the vicinity of the 
transition region between the plateaus, when breaks 
down, the time-dependence mixes states with N — no 
and N — no + l. This means that the width of the transi- 
tion region is given by Te/ / , and in the zero-temperature 
limit saturates to 1/t. In this regime the nearly adia- 
batic picture fails. The width of the charge states due to 
tunneling T{t) decreases with time. When T{t) ^ l/r, 
the system can no longer follow the changing tunneling 
coupling. Effectively, it can be described within the sud- 
den approximation, where T{t) is replaced by the step- 
function, r{t) rsO{—t), Ts ^ 1/t . The occupation of 
the dot at t ^ CO is therefore determined by (|[^) with 
time-independent tunneling amplitude, corresponding to 
the width Ts- 

Due to the interaction term in (H), the model (||-^ is 
still difficult to solve analytically. To simplify the dis- 
cussion, we limit our attention to the interval of J\fg, 
which includes only one transition region between the 
plateaus: no < Mg < no + 1. Furthermore, we neglect 
the spin degeneracy, and consider the limit of the large 
single-particle level spacing in the dot, A 3> l/r. With 
these restrictions, at low temperature T ^ A, only the 



lowest energy configurations, corresponding to A^ = no 
and A^ = no -l- 1, are important. Since these states are 
non-degenerate, one can introduce the fermion operator 
d = I no) (no -|- 1| to describe transitions between these 
states, and replace (0) by 



Hdot = EoSd, Eo = 2E, (1/2 + no - Afg) . 



(9) 



The advantage of the model (||),(H),(|^) is that it is exactly 
solvable for arbitrary V{t). Indeed, the occupation of the 
dot n (t) = i^d'^ {t)d{t)') = N (t) — no satisfies equation of 
motionB 



dt 



n{t) = ~T{t)n{t) + / denF{e)A (e, t) , 



(10) 



A{e,t) = Im / dt' ^T{t)T{t')e'^^'-*'^G"' {t,t') . (11) 



Here T{t) = 2ttpV^ (t) is the width of the charge state, 
z/ is density of states of conduction electrons at the 
Fermi level, np{e) is Fermi function and G^{t,t') = 
— 10 (t — t') ({(i(i), (i^(t')}) is the exact retarded Green 
function of the dot: 

G^- {t, t') = -ie {t - t') e-' '^*i[^o-»r(*i)/2] 



Solution for A'o follows by simple integration of (|lO|-[llD 
With V [t) given by Eq. (0), the result is 



Ao - no 



To_ 
27r 



de- 



npis) 



(12) 



cosh [(e — Eo) To] ' 
where To = ttt /2. At zero temperature, (|l2|) reduces to 



A^n 



no = -tan^i [e'^''^°] 



(13) 



At finite temperature, (12) is described very well by the 
Fermi function 



A^o - "0 
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(14) 



with an effective temperature Tg// = v (c/to) -tT^- We 



found that c = 0.88 gives the best numerical fit) 

For {Afg — > no), Eqs. (|l|) and (|lj) give the following 
expression for the slope of the plateau: 



h \dN, 



{2E,/T,ff)e 



(15) 



Here /q = noef corresponds to perfect quantization. 
Strictly speaking, to obtain the correct value of the slope 
precisely in the middle of the plateau, the two-state ap- 
proximation is not sufficient: state with A^ = no — 1 makes 
exactly the same contribution, as that with A^ = no -I- 1. 
This complication, however, should not affect signifi- 
cantly the validity of (|T5|): the exact result differs from 
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( p^ ) by the factor of the order of 1 onlyM. Due to the ex- 
ponential factor in (|l5|), S depends very strongly on the 
ratio Ec/Teff. For example, for Ec/Tefj = 10, S* - 10"^^ 
while for EjTeff = 20, S ^ 10"'^. 

According to the discussion above, in the transition 
region, the result can be obtained with the sudden ap- 
proximation. Thus, we have 

AT [a ( r,/27r 

A'o - no « / denple) ^ , 

or, for T = 0, TVo - no « 1/2 - (l/vr) tan~i {2Eo/rs). 
This expression indeed coincides with (|l3|) in the limit 
Eq 0, if Ts/2 = ro~^ Note that the width of the tran- 
sition region is determined by essentially the same Tg/ / , 
that enters ( [15|). 

The model (Hjq) introduced above allows study of the 
nonadiabatic effects at the short time-scale. As the sys- 
tem evolves with time, the SAW-induced potential mini- 
mum moves uphill (see Fig. 2), and may eventually cross 
the Fermi level. Due to the residual tunneling coupling 
in this regime, the leakage from the dot will introduce ad- 
ditional corrections to (|^). These corrections, however, 
do not affect strongly the slope of the plateaus , and 
can be taken into account by multiplying by the leak- 
age factor Pi < 1: I = PiNoef, where Pi is expected to 
depend on system parameters, such as gate voltage and 
SAW power. Thus, the exact value of the quantized cur- 
rent Iq — PiUgef does not necessarily coincides with the 
transfer of precisely integer number of electrons per pe- 
riod no, and the plateaus can move in parameter space. 

In conclusion, we have shown that at low tempera- 
ture, long-range Coulomb interactions may have dra- 
matic effect on the accuracy of the quantization of 
the single-electron transport driven by surface acoustic 
waves through a narrow constriction, formed in two- 
dimensional electron gas. The effect of screening on 
the SAW-induced potential near the edge of 2DEG can 
be described by a single parameter - the time r of the 
switching-off of the tunneling coupling between 2DEG 
and the moving quantum dot. As a result, both the slope 
of the plateaus and the width of the transition regions 
between the plateaus saturate at low temperature to the 
values, determined by the characteristic energy scale for 
nonadiabatic corrections h/r. 
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